Nonperiodic Flux to Voltage Conversion of an Arithmetic Series Array of dc SQUIDs 
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A theoretical study on the voltage response function (V) of a series array of dc SQUIDs is presented 
in which the elementary dc SQUID loops vary in size and, possibly, in orientation. Such series 
arrays of two-junction SQUIDs possess voltage response functions vs. external magnetic field B 
that differ substantially from those of corresponding regular series arrays with identical loop-areas, 
while maintaining a large voltage swing as well as a low noise level. Applications include the design 
of current amplifiers and quantum interference filters. 



Recently series arrays of dc SQUIDs have been suc- 
cessfully exploited as current amplifiers with wide.band- 
width, large dynamic range and low noise levelEJ. Us- 
ing thin-film Nb-technology, amplifiers consisting out of 
up to 10 3 identical dc SQUID loops have been fabri- 
cated. Such serial devices are characterized by large volt- 
age swings of several mV and current-to- voltage transfer 
functions of some V/ mA so that a direct connection to a 
room temperature preamplifier is feasible. 

At present only series arrays consisting out of identical 
dc SQUID loops have been described in the literatureEl. 
The dc voltage response function (V) of such regular se- 
ries arrays displays a <I>o-periodicity just as a single dc 
SQUID where $o = is the elementary flux quantumo. 
However, this periodicity may represent a serious limi- 
tation for the modes of operation of such devices. To 
get a linear flux-to-voltage conversion a feedback circuit 
is used in most applications in which the dc SQUID 
acts as a null-flux detector. Furthermore special elec- 
tronic devices and efficient background shielding is of- 
ten requireda. Here theoretical studies on the voltage 
response of arithmetic and irregular series arrays of dc 
SQUID loops are presented, in which the individual loop- 
areas are not all equal. For such devices the advantages 
of dc SQUID series arrays are preserved but limitations 
due to the periodicity of the voltage response can be cir- 
cumvented. 

The arrays under consideration consist of N two junc- 
tion SQUID loops connected in series. The bias current 
h is fed into the array as indicated in Fig. (|l|) . In general 
the areas of the N loops in a generic array differ in size 
and, possibly, in orientation. Let a„ be the orientated 
area element of the n th dc SQUID loop. The magnetic 
flux threading a n is then $ n = (B, a„), where B is the 
total magnetic field. Taking into account all inductive 
effects in the dc SQUID array the total magnetic field 
B = B^ + B^ s ) + B^ c ) is a superposition of the primary 
magnetic field B^ p ' one wants to measure, the secondary 
magnetic field B^ s ^ induced by currents that flow in the 
dc SQUID array and, possibly, the magnetic compensa- 
tion field B( c ) induced by compensation current(s) I ramp 
flowing through a set of suitable orientated compensation 
coils or wires. 

As a special case a generic planar dc SQUID array is 
shown in Fig.(g). The primary magnetic flux is coupled 
into the individual dc SQUID loops by the signal or in- 



put current Ii„ p which flows through a common input 
flux bias line. In typical applications Ii np is the current 
provided by a pick-up loop or by other signal sources. 
In addition the flux in the individual dc SQUID loops 
may be modulated by the current I CO mp flowing through 
a common compensation or control flux bias line. 
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FIG. 1. Schematic diagram of a generic planar series array 
of two-junction SQUIDs. a n denotes the orientated area ele- 
ment of the n th dc SQUID loop. The Josephson junctions are 
indicated by crosses (x) . The array is biased with the current 
h and produces an dc output voltage (V) . The magnetic flux 
is coupled into the dc SQUID loops by Ii np flowing in a com- 
mon input flux bias line and also by I CO mp in a control flux 
bias line. 



Let Mn an d Mn 1 denote mutual inductances of the 
currents I CO mp and Ii np associated with the n th loop with 
area-element a n and let be the inductance matrix de- 
scribing all other inductive couplings in the circuit!! The 
total magnetic flux in the n th loop is then given by <f>„ = 
M n p) I inp + M n c) I comp + L A I, where I = {h , • ■ ■ , h n) is 
the vector of the 2N currents I n through the Josephson 
junctions. Since I CO mp should be able to compensate the 

(p) 

primary induced magnetic fluxes M n Imp in all loops 

(c) (v) 

simultaneously the relation M n = a M n should hold 
for all n, where a is a constant. _ 

Within the range of validity of the RCSJ-modelQ the 
current I m through the m th array junction with gauge in- 
variant phase difference ip m , critical current J c , m , normal 
resistance R m , and capacitance C m is given by 
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For simplicity identical parameters for the junctions are 
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supposed in the following discussion: I c ,m — la Rm = 
R and C' m = C. The dimensionless parameter (3c = 
27r/ |> fl ° characterizes then the capacitance effect of the 
junctions within this modelR 

In the overdamped regime, C = 0, and for vanishing 
inductive coupling in the array, ~La — 0, an analytical ex- 
pression for the dc voltage response function (V) of the 
dc SQUID series can be determined. Under conditions 
where a constant current lb > 2 I c is biased, and assum- 
ing for simplicity a static magnetic field B = B' p ^ +B' c \ 
the output voltage (V) n of the n th dc SQUID in the ar- 
ray as a function of^hc bias current lb and of the applied 
flux <!>„ is given by 
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The voltage response (V) of a series array containing N 
dc SQUIDs is then given as the sum of the voltages (V) : 



(3) 



Since the maximum voltage swing increases linearly with 
the number N of dc SQUIDs in the array the gain §j^- 
is proportional to N. However, the noise in the array 
increases only as , because the noise voltages implied 
by the junctions in the individual SQJJID-loops of the 
array are expected to add incoherentlyu. These features 
make generic series arrays attractive for amplifiers and 
magnetometers . 

Consider, as a special case, a regular array, consisting 
out of N identical dc SQUID loops so that a n = ajv with 

identical mutual inductances M, 

(p) 



M^ c) and M 
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N of the flux bias lines. As a result the flux threading 



M 

each individual loop is then given by <!>„ = $n. In this 
case the voltage response, (V) = N (V) N , displays as a 
function of flux a <E>o-periodicity. 

A voltage response function (V) with a much longer 
period may be obtained in an arithmetic dc SQUID series 
array where the orientated area elements increase in size 
according to the arithmetic relation 



a„ = — ajv, ra=l, 
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Provided the mutual inductances are chosen such that 



M 



and 



the fluxes are 



|a„| 

|ajv| -- J jV lv± n — | ajv | lv± N 

determined by $„ = ■§> $jv, where is the flux through 
the largest area element a^v • In this case the individual dc 
SQUIDs operate with different periodicity in flux. This 
implies that the voltage response (V) of the whole arith- 
metic array as a function of $ at is periodic with period 
N$o. This period is N times larger than for a regular 
array, consisting of N identical dc SQUID loops with size 
a N . 



In Fig.fg) the voltage response for an arithmetic se- 
ries array with N — 100 dc SQUIDs is shown. The 
bias current lb was adjusted for maximum voltage swing 
slightly above 2I C . In the depicted interval the voltage 
response shows pronounced antipeaks at $n = and at 
<I>Ar = N $o- Between these two peaks the output voltage 
oscillates with a very small amplitude around an average 
voltage drop V. 
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FIG. 2. Voltage response (V) according to Eq.(g) in units 
of I c Rvs. external flux $jv through the largest area element 
ajv for an arithmetic series array with N = 100 dc SQUIDs 
for bias current /{, = 2.2 7 C and vanishing inductive coupling 
{PL,max = 0). The dc SQUID loops increase in size according 
to the linear relation a„ = ajv- In the inset the theoreti- 
cal curve (solid line, Eq.@) is compared with the analytical 
approximation (circles, Eq.(^)). 

For a sufficiently high number of elements in the arith- 
metic array the flux-distribution over the dc SQUID 
loops is uniform between some $> m in and <fr m ax- Let 
k = ® min be the flux spread coefficient, / = n $ g° x 
the frustration of the greatest loop in the array and let 
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h 

2 1, 



1 be the characteristic frequency of a 



single dc SQUID for zero flux. If the frustration / is 
limited to the periodicity interval |/| < 7r ^ the voltage 
response (V) as a function of / can be approximated by 
an integral: 
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where m = — u>q and Elf, m) denotes the elliptic in- 
tegral of the second kindj. The voltage versus flux re- 
lation displayed in Fig.(||) shows an excellent agreement 
between the theoretical result Eq.(||) and its analytical 
approximation Eq.(||). For bias currents lb = 2.2 I c there 
results an average voltage drop of V = 0.81 I C RN. 

At first sight all the above results apply only in the 
overdamped junction regime and only for vanishing in- 
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ductive couplings in the arithmetic dc SQUID array, i.e. 
for Pc = and Lia = 0. A dimensionless measure for the 
inductive effects in a single dc SQUID loop with loop- 
inductance L is (3l — ^gk. Dependent on /3j, there exists 
an optimal size for any dc SQUID loopH which should co- 
incide with the maximal loop size a. max in a series array. 
The corresponding maximal Ph,max is then a measure of 
the self- and mutual-inductive couplings among the cur- 
rents flowing in the dc SQUID array. 

Taking into account all inductive couplings in 
the arithmetic array our computer simulations of the 
full non-linear dynamicsa of the 2 TV coupled Josephson 
phases tpi reveal that expression Eq.(Q) compares quali- 
tatively well with the calculated voltage response func- 
tions. In particular the long periodicity of the voltage 
(V) vs. flux <&at relation is not affected by inductive 
effects. Qualitatively similar behavior was found for a 
parallel multi-junction interferometenl 

As far as the irregularity is concerned, we also find 
that in an arithmetic array (V) is very responsive to 
adding small random fluctuations to the size distributions 
of the area elements. In this case (V) becomes nonperi- 
odic with a pronounced antipeak only around $jv = 0. 
For fixed Pc, but different strengths of the magnetic cou- 
pling 0L,max, Fig-® shows the voltage response for an 
arithmetic planar series array with N = 100 dc SQUIDs 
where small random fluctuations were added to the loop 
size distribution Eq.(||). 

For fixed Pc the global minimum of (V) at <E>jv = 
depends only slightly on the strength Ph,max of the mag- 
netic coupling whereas the voltage branch (V) « V de- 
pends strongly on Pr,,max, as can be seen in the inset 
of Fig.([|). For increasing Ph^max the voltage swing and 
the maximum of the voltage transfer function »^ de- 
crease as it is the case for conventional single dc SQUIDs 
when (3l increasesa. If on the other hand PL.max is kept 
constant the global minimum of (V) at $jy = is an 
increasing function of Pc as a comparison of Fig. (^) and 
Fig.(||) reveals. In addition Fig.(|p,) shows an enhance- 
ment of the branch (V) w V for Pc = 0.5 in comparison 
to the underdamped case, since for single dc SQUIDs 
the effect of the capacitance pc is to increase the volt- 
age near integer values of the applied fluxcl. Although 
the total area A tot = |ajy| of an arithmetic array 

is smaller by a factor of 2 compared to a regular array 
with A to t — N |ajv|, this effect provides that the volt- 
age transfer function of both arrays is comparable. This 
suggests a higher integration density for arithmetic array 
circuits on chip. 

In summary, arithmetic series arrays of two-junction 
SQUIDs possess voltage response functions with a much 
longer period with respect to the applied magnetic field 
than regular series arrays, while maintaining a compara- 
ble transfer function and a low noise level. In particular 
if the loop-sizes or loop-orientations are distributed ran- 
domly the voltage response becomes nonperiodic with 
a single pronounced antipeak only around zero applied 



magnetic field. These features are preserved when all 
capacitive and inductive effects are taken into account. 
Therefore arithmetic or irregular series arrays of dc 
SQUIDs can be used as quantum interference filters for 
various applications, including, e.g., the relatively sim- 
ple and extremely sensitive measurement of the absolute 
strength of magnetic fields. 
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FIG. 3. Voltage response (V) in units of I c R vs. exter- 
nal flux <i>jv for an arithmetic series array with N = 100 
for It — 2.2 I c , Pc = 0.5, and Phmax = 0.3. Small ran- 
dom fluctuations were added to the loop size distribution 
an = j? sln- The inset depicts (V) for various inductive cou- 
plings around the global minimum at "1>jv = 0: a) PL,max = 0, 
b) PL,max = 0.3, and c) PL, ma x = 0.7. 
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